Critical Exponents for the Metal-Insulator 
Transition in Disordered Systems 



Angus MacKinnon 
Blackett Lab., Imperial College, London SW7 2BZ 

February 1, 2008 



Short title: Metal-Insulator Transition 
PACS number(s): 71.30.+h, 71.55 Jv, 72.15Rn 
Submitted to: Journal of Physics: Condensed Matter 

Abstract 

The critical exponents of the metal-insulator transition in disor- 
dered systems have been the subject of much published work contain- 
ing often contradictory results. Values ranging between 1/2 and 2 can 
be found even in the recent literature. In this paper the results of a 
long term study of the transition are presented. The data have been 
calculated with sufficient accuracy (0.2%) that the calculated expo- 
nent can be quoted as s = v = 1.54 ± 0.08 with confidence. The 
reasons for the previous scatter of results is discussed. 



1 Introduction 



The metal-insulator transition in disordered systems has been the subject 
of theoretical and experimental work at least since Anderson (|i~958l) . The 
similarities with thermodynamic phase transitions had been noted by several 
authors ( Thouless, 1974j ; [Wegner, 1976| ) but it was not until 1979 that a 
usable formulation of the renormalisation group or scaling theory became 
available ([Abrahams et al, 1979 ; VVcgncr. 1979 ; [Efctov, 1983 ). The basic 
assumption of these theories, that the behaviour could be described by a 
single parameter scaling theory, was confirmed in numerical calculations by 
the present author ( |MacKinnon fc Kramer, 1981| ; |MacKinnon fc Kramer,| 
1983| ) . For a recent review of the area see Kramer and MacKinnon ( |1994| ) . 



In spite of the progress made the exponents, s and u, describing the 
behaviour of the conductivity and the localisation length respectively have 
proven difficult to calculate reliably. For some time there appeared to be a 
consensus between theory and experiment that both exponents were equal 
to unity, but more recently this has been called into question from both the 
theoretical (e.g. ( [Kravtsov fc Lerner, 198"^ ; [Lerner, 1991| ) ) and from the 
experimental ( (Stupp et al, 1993|) side. 

Numerical results have been scattered at least between 0.5 and 2 with 
numerous attempts at developing alternative methods of calculation. A good 
example of the difficulties is given by the contrast between calculations for 



the Anderson model with rectangular or Gaussian disorder (|Kramer et al. 



1990|) . Using identical methods the exponents obtained were about 1.5 and 



1.0 for the rectangular and Gaussian distributions respectively. It is clearly 
unreasonable for the exponents for these two cases to be different. In fact if 
they were different then it would call into question the justification of the use 
of any simple model Hamiltonian to describe the transition and so undermine 
the whole foundation of the subject. 

In this paper the results of calculations carried out over several years are 
presented. All the basic results have an accuracy of at least 0.2% which 
enables the critical exponents to be calculated much more accurately than 
when the conventional 1% is used. 
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2 Transfer Matrix Calculations 



The transfer matrix method has been discussed in numerous papers ([MacK- 



innon fc Kramer, 1983| ; |Pichard fc Sarma, 1981| ) so only the briefest outline 



will be attempted here. 

The starting point is the usual Anderson ( |1958|) Hamiltonian 

H = ^2e i \i><i\+Y i V ij \i>< j\ (1) 

where = Vq between nearest neighbours on a simple cubic lattice and zero 
otherwise. In this work Vq = 1 is chosen and will therefore not be mentioned 
explicitly. The diagonal elements e« are independent random numbers chosen 
either from a uniform rectangular distribution with —1/2W < £j < + 1/2W 
or from a Gaussian distribution of standard deviation a. For purposes of 
comparison between the two cases an effective W for the Gaussian case may 
be defined by equating the variances as W 2 = 12a 2 . 

In terms of the coefficients a, of the wavefunctions on each site the 
Schrodinger equation may be written in the form 

Ecii = ejdj + ^a,-. (2) 

Consider now a long bar composed of L slices of cross-section M x M. 
By combining the a,s from each slice into a vector Aj (§) can be written in 
the concise form 

EA n = H n A n + A n+1 + A n _i (3) 

where the subscripts n now refer to slices and matrix H ra is the Hamiltonian 
for slice n. By rearranging (^) the transfer matrix is obtained 

(4a) 
(4b) 
(4c) 
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A theorem attributed to Oseledec ( |1968| ) states that 

lim (T|,T n ) 1/n = M (5) 

where M is a well defined matrix and T n are products of random matrices. 
The logarithms of the eigenvalues of M are referred to as Lyapunov exponents 
and occur in pairs which are reciprocals of one another. By comparison with 
(|J) the Lyapunov exponents may be identified with the rate of exponential 
rise (or fall) of the wave functions. In fact the smallest exponent corresponds 
to the longest decay length and hence to the localisation length of the system. 

In principle then it is necessary to calculate T n for large n, and diagonalise 
T'T. Unfortunately the calculation is not quite so simple: the different 
eigenvalues of T^T rise at different rates so that the smallest, which we 
seek, rapidly becomes insignificant compared to the largest and is lost in the 
numerical rounding error. Typically this happens after about 10 steps. 

2.1 Orthogonalisation 

In order to obtain the smallest Lyapunov exponent it is necessary to overcome 
this loss of numerical significance. This can be achieved in more than one 
way of which the orthogonalisation method is employed here. 

After about 10 matrices have been multiplied together the columns of 
the product matrix are orthogonalised to each other and normalised. This is 
equivalent to multiplying the product from the right by an appropriate ma- 
trix. This orthonormalisation process automatically separates the different 
exponentially growing contributions. 

The process is repeated every 10 or so steps and the logarithm of the 
length of the vector closest to unity is stored. The Lyapunov exponent is 
given by the mean value of these logarithms divided by the number of steps 
between orthonormalisations. In practice it is necessary to use only 50% or 
M x M vectors rather than the full 2 x M x M as the required vector is 
invariably the M x Mth. 

The error in the Lyapunov exponent can be estimated from the variance 
corresponding to the mean exponent. Although this estimate could be biased 
by correlations between the different contributions this is not found to be a 
serious problem in practice, at least when the localisation length is short 
compared with the distance between orthogonalisation steps. 
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The optimum frequency of orthogonalisation steps can be estimated by 
comparing the length of the M x Mth vector before and after orthogonalisa- 
tion. The ratio should not be allowed to get close to the machine accuracy. 



3 Scaling Theory 

The inverse of the smallest Lyapunov exponent is the localisation length 
Am- The renormalised length A = Xm/M is found to obey a scaling theory 
([MacKinnon fc Kramer, 1981| ; [MacKinnon fc Kramer, 1983| ) such that 



din A „ , s . . 

which has solutions of the form 

A = f (M/0 (7) 

where £ is a characteristic length scale which can be identified with the 
localisation length of the insulator and which scales as the reciprocal of the 
resistivity of the metallic phase ( [MacKinnon fc Kramer, 1983|) . 



In 3D (|) always has a fixed point x = which corresponds to the metal- 
insulator transition. The behaviour close to the transition can be found by 
linearising (||) and solving to obtain 

In A = In A c + A(r - r c )M a (8) 

where r is the disorder W or a, A c and r c represent the critical A and 
disorder respectively, and A and a are constants. By comparing (|^) and 
(HD an expression for £ can be obtained in the form 

i-\r~r c \ 1/a (9) 

so that the localisation length exponent v is given by v = 1/a. Since it 
is well known ( |Wegner, 1976| ; [Abrahams et al, 1979| ) that the conductivity 
exponent s is related to v by s = [d — 1)v then by fitting (|8]) to the data and 
calculating a both exponents can be obtained. 
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3.1 Deviations from Scaling 



One simple feature of is that, when In A is plotted against r, the curves 
for different M intersect at a common point (lnA c ,r c ). In practice the data 
do not behave in exactly this way. There is a small deviation from scaling. 
This deviation could be taken into account by adding an extra term to (H) 
which depends on M but not on r. Consider, however, the form 

In A = ArM a + B (M) (10) 

which represents the most general form of such a correction. If a specific 
form for the correction were assumed it would require at least 4 independent 
fitting parameters to represent B(M), including A c and r c , and may still not 
represent the true deviation from scaling. It seems better therefore to fit an 
independent B(M) for each value of M and therefore to make no assumption 
about the nature of the deviation from scaling, other than that it is non- 
critical, and therefore independent of r, in the region of interest. By fitting 
the data to (P~D| ) in this way the exponent a is derived solely from the gradient 
of In A vs. t and the intercept is allowed to float. The results of such fits are 
shown in figure [l]. 



3.2 Data Fitting 

The data can be fitted to ([U]) by iteratively using a standard least squares 
procedure. Care is required with the non-linear parameter a. The quality 
of the fit can be tested by computing \ 2 defined as 

x 2 = y^ {MM? + B(M i )-lnA i ) 2 (n) 

where % runs over all data points and cr, is the error in point i. After fitting 
X 2 should be approximately equal to the number of data points less the 
number of fitted parameters. Hence the value of \ 2 provides a measure of 
the quality of the fit. In the results presented here the range of values of 
disorder round the critical value was chosen such that y 2 conforms to this 
condition. Then a large number of additional points was calculated inside 
this range. An important side effect of this procedure is that the apparently 
acceptable range of disorder around the fixed point gets narrower as the 
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calculations become more accurate. It is therefore important to test whether 
any apparent change in the fitted exponent is due to this narrowing. 

The values of the ideal and the fitted x 2 as we U as the range considered 
are shown in table [|. Using 4 < M < 12 and the widest range of disorder 
s = v = 1.53 ± 0.04 and s = v = 1.48 ± 0.05 for rectangular and Gaussian 
cases respectively. 

3.3 Statistical and Systematic Errors 

The statistical error in the fitted critical exponent is easily estimated from 
the least squares fitting procedure. Systematic errors are more difficult to 
take into account. In this work an attempt is made to consider 3 sources of 
systematic error: 

• Limited range of system sizes: 4 < M < 12 has been considered and 
the effect of ignoring the smaller system sizes tested. 

• Width of the critical region: the maximum range of disorder is imposed 
by x 2 but may still be too large. The effect of narrowing this range 
still further has been tested. 

• The choice of distribution of random numbers: this has been tested by 
comparing the rectangular and Gaussian cases. 

These tests are represented in figure |^. Unfortunately the general increase in 
the error bars due to ignoring data tends to mask any systematic changes. 
There does however appear to be a general increase in the exponents when the 
M = 4 data is eliminated and a tendency for the Gaussian data to lie below 
the rectangular. From this data s — v ~ 1.54 ± 0.08 has been estimated, 
where the error bar may be somewhat wider than necessary. 

4 Results and Conclusions 

The results are summarised in table |l|. All these results have been calculated 
in the middle of the band (i.e. E = 0), but there is ample evidence that for the 
models considered here, this point is not special and is truly representative 
of the whole band, at least in the range —6 < E < 6. 
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Rectangular 


Gaussian 


Exponent 
Disorder Range 
System Sizes 
X 2 (expected) 
X 2 (fitted) 

w c 

A c 


1.515 ±0.033 
16.2 < W < 16.8 
4 < M < 12 

142 

126 
16.50 ±0.05 
4.763 ±0.015 
0.580 ±0.005 


1.484 ± 0.048 
21.0 < W < 21.5 
4 < M < 12 
97 
75 

21.20 ±0.06 
6.120 ±0.018 
0.580 ±0.005 



Table 1: N.B: The estimates of W c and A c are based on the values given by 
several different fitting procedures. 



Unlike previous calculations ( [Kramer et al, 1990| ) the exponents calcu- 
lated for the two distributions now overlap well and are therefore consistent 
with the common assumption that simply changing the distribution does not 
change the universality class and hence the critical exponent. The discrep- 
ancy reported previously is presumably due to insufficient accuracy in the 
raw data and consequent assumption of a critical range of disorder which was 
too wide. 

This may have consequences for experiment as it seems to suggest that 
it is possible to obtain an exponent of unity simply by using too wide a 
range of data around the critical disorder, energy, pressure, etc. It should 
also be borne in mind that the influence of interactions may also account for 
differences between experimental results and those based on a model of non- 
interacting electrons. For this reason it may be more realistic to compare the 
present results with photonic or acoustic rather than electronic experiments. 

In summary, the critical exponent of the Anderson model of the metal- 
insulator transition is s = v = 1.54 ± 0.08. 

Acknowledgements 

This work has profited from many useful discussions with B.Kramer, M.Schreiber, 
J.B.Pendry, P.M.Bell, R.B.S.Oakeshott, E.A.Johnson, and P.J.Roberts. The 
financial support of the UK SERC and the European Union, through SCI- 
ENCE grant SCC*-CT90-0020, is gratefully acknowledged. 



7 



Figure Captions 

1. A vs. W, for (a) rectangular and (b) Gaussian distributions. The data 
are represented by dots with differing symbols for different system sizes 
with 4 < M < 12 increasing in the direction of the arrow. Each point 
is accurate to 0.2%. The lines are fitted using fllPf) . 

2. Fitted critical exponents for rectangular (Diamonds) and Gaussian 
(Squares) distributions. The abscissa? represent the smallest system 
size taken into account (with small offsets for clarity). In each group 
the width of the fitted region is (from left to right) (16.2 < W < 
16.8) -> (16.3 < W < 16.7) -> (16.4 < W < 16.6) and (21.0 < W < 
21.5) -> (21.05 < W < 21.45) -> (21.1 < W < 21.4) for rectangular 
and Gaussian cases respectively. The dotted lines represent the range 
s = v = 1.54 ±0.8. 
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